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Abstract 

The previously introduced Quantum Arnold Transformation, a unitary operator mapping 
the solutions of the Schrodinger equation for time dependent quadratic Hamiltonians into 
the solutions for the free particle, is revised and some interesting extensions are introduced, 
providing in particular a generalization of the Ermakov-Pinney equation. 

1 Introduction 

In [TJ, V.I. Arnold introduced a theorem stating that the family of graphs of solutions of any 
linear second-order differential equation (LSODE), with arbitrary time dependent coefficients, 
is locally diffeomorphic to the family of graphs of solutions of the simplest one-dimensional 
equation of motion, i.e., the equation of motion of a free particle. This important result allows 
to establish a direct relationship between the symmetries of an arbitrary classical linear system, 
whose equation of motion is a LSODE, with that of the free particle, which is known to have 
the largest possible symmetry group, SL(3,M). The diffeomorphism was explicitly given in [1], 
namely the Classical Arnold Transformation (CAT), and can be seen as a particular class of 
more general Lie transformations relating systems defined by second order differential equations 
with the free particle [2]. 

In [3], the authors, in collaboration with V. Aldaya and F. Cossfo, introduced the quantum 
version of the CAT, the Quantum Arnold Transformation (henceforth QAT), as a unitary map 
that relates the Hilbert space of solutions of the time-dependent Schrodinger equation for a 
Generalized Caldirola-Kanai oscillator (the quantum version of a classical system whose equation 
of motion is a LSODE) into the corresponding one for the free particle. In the present paper 
we revise the QAT and extend it in two respects. First, we give the explicit transformation 
connecting quantum physical systems with arbitrary time-dependent quadratic Hamiltonians 
with the free particle system, here called the Gauged QAT. Second and more importantly, we 
describe the transformation that links any LSODE-system to any other LSODE-system in terms 
of CAT's, giving rise to a generalization of the Ermakov-Pinney equation, and for that reason 
we shall call it the Arnold-Ermakov-Pinney transformation, providing also its quantum version. 

Although the QAT is very recent, the CAT is more than thirty years old and, since it was 
introduced, many authors have used it (even without realizing it) and have built something 
similar to the QAT. Even before that, Lewis and Riesenfeld (1969) [1] introduced a technique 
to obtain solutions of the time-dependent Schrodinger equation (TDSE) for a time-dependent 



quadratic Hamiltonian (TDQH) as eigenfunctions of quadratic invariants. For that purpose 
they wrote the solutions in terms of auxiliary variables that satisfy the classical equations 
of motion (something that resembles the CAT). Dodonov & Man'ko (1979) [5] constructed 
invariant operators for the damped harmonic oscillator and introduced coherent states, using a 
method similar to that of Lewis and Riesenfeld. Jackiw (1980) |6] gave (implicitly) the quantum 
transformation from the harmonic oscillator (even with a 1/x 2 term) to the free particle when 
studying the symmetries of the magnetic monopole. Junker & Inomata (1985) |7J gave the 
transformation of the propagator, in a path integral approach, for a TDQH, into the free one 
(the equivalent of the QAT, but in terms of propagators). Takagi (1990) [8] gave the quantum 
transformation from the harmonic oscillator to the free particle, interpreted as the change to 
comoving coordinates. Bluman & Shtelen (1996) |9j gave the (non-local) transformation of 
the TDSE for a TDQH plus a non-linear term into the free particle one, in the context of 
transformations of PDEs. Kagan et al. [10] and independently Castin & Dum |llj (1996) 
introduced a scaling transformation in the Gross-Pitaevskii equation describing Bose-Einstein 
Condensates (BEC) which is related to the QAT. Suslov et al. (2010) [12] computed the 
propagator for a TDQH using the classical equations. 

To the best of the authors' knowledge, the Quantum Arnold Transformation can be not 
only very useful to perform certain analytic calculations, but key to unify different concepts 
scattered through the mentioned literature in a simplified framework. 

The paper is organized as follows. First, in Section [2] we revise the Classical Arnold Trans- 
formation. After that we go on with the Quantum Arnold Transformation in Section [3j Finally 
two extensions are described in Section [4j the Gauged QAT and the Arnold- Ermakov-Pinney 
transformation. 

2 The Classical Arnold Transformation 

To set up the framework in which the classical Arnold transformation acts, let us briefly remind 
the subject of Lie symmetries of ordinary differential equations. 

2.1 Lie transformations 

A Lie symmetry of an ordinary differential equation (ODE) is a coordinate transformation that 
sends solutions into solutions. The problem of determining the Lie symmetries of an ODE is 
rather old, and S. Lie gave the main results at the end of the nineteenth century [2]. One 
of these results was that a second order differential equation (SODE) y" = F(x, y, y') has the 
maximal number of Lie symmetries (SX(3,M)) if it can be transformed to the free equation by 
a point transformation: 

V =F{x,y,y) =>• y =0. (1) 
This linearization is possible if the ODE is of the form: 

y" = E 3 {x,y)(y'f + E 2 (x,y)(y') 2 + E 1 (x,y)y' + E {x,y) , (2) 

with Ei(x,y) satisfying some integrability conditions (see, for instance, [T3"l [H]). 

There is a nice geometric interpretation of this condition in terms of projective geometry. 
The non-linear SODE ^ is obtained by projection from the geodesic equations in a 2-dim 
Riemannian manifold. The coefficients Ei(x,y) are in one-to-one correspondence with Thomas 
projective parameters n, and the integrability conditions that they satisfy are the conditions 
for the Riemann tensor to be zero (see [1'6\ I14|). 
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V.I. Arnold named this process rectification or straightening of the trajectories, and studied 
the case of linear SODE (LSODE), giving explicitly the point transformation for this case pp. 
In the next subsection we describe it in detail. 



2.2 Classical Arnold transformation 

A General Linear Second Order Differential Equation (LSODE) is given by the differential 
equation: 

x + fx + u?x = A, (3) 

where f,u and A are functions of t. The Classical Arnold Transformation (CAT) is a point 
transformation that is a local diffeomorphism: 

A: IxT ^IxT / r = g| = Jj 

where T and T are, in general, open intervals containing to and 0, respectively, u\ and U2 are 
independent solutions of the homogeneous LSODE satisfying the canonicity conditions: 

«i(*o) = = u' 2 (to) , u' x (to) = 1 = u 2 (to) , (5) 

Up is a particular solution of the inhomogeneous LSODE satisfying u p (to) = u'(to) = 0, and 
W(t) = u\U2 — u\u 2 = e~f is the Wronskian of the two solutions. Here to is an arbitrary time, 
conveniently chosen to be to = (see [3] for details). 

The CAT transforms the original LSODE ^ into that of the free particle, up to a factor: 

o A W 

x + fx + oj 2 x = A — ► -=■ R = 0. (6) 

The presence of this factor implies that patches of trajectories of ^ are transformed into 
patches of straight (free) trajectories. In fact, an arbitrary trajectory solution of ^ can be 
written as x(t) = Au\(t) + Bu 2 (t) + u p (t), and the CAT sends it to k(t) = At + B. While t 
varies in the interval T defined by two consecutive zeros of u 2 (t) (containing to), r varies in the 
range of the map defined by ^ft) • ^ n ^ e case ^ n "^^icl 1 u 2(t) has one zero, T is (left- or right-) 
unbounded, and, if it has no zeros, T is M. 

Even though the CAT is a local diffeomorphism, it can be defined for an arbitrary time 
to- Thus different CATs can be defined for different times to and cover in this way a complete 
trajectory of We shall show with the example of the harmonic oscillator how this can be 
done. 



2.2.1 The example of the harmonic oscillator 

The harmonic oscillator (HO) is the best example to understand how the CAT works. For this 
case, and considering A = 0, the two solutions are u\(t) = ^ sin(wt) and u 2 (t) = cos(wt). The 
open interval T defined by two consecutive zeros of u 2 (t), and containing to = 0, is (— #j), 
and the CAT A and its inverse A -1 are then written as: 

r =^i = ^tan(u,i) (7) 

(8) 




U2(i) cos(arf) ' U2(t) 

x = cos(arctan(wr))K , t = ^ arctan(wt) . 

K 

~ Vl+L0 2 T 2 
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In this case t G R. Pictorially, the CAT for the HO can be represented as in Figure [TJ 
where velocities have also been included in the graphic for clarity. Here A maps the solid 
part of the helix (half a period of a harmonic oscillator trajectory) into the whole line (a free 
particle trajectory). The horizontal plane represents the space of all possible initial conditions 
at t = = t. Note that both trajectories are tangent when projected onto this plane, due to 
the conditions (pi). See [TH] for more details in this case. For the CAT to map other patches of 



/ 



(k(t),k(t),t) 




Figure 1: Depiction of the CAT for the harmonic oscillator (adapted from [15j). 

the HO trajectories into the free particle trajectories, different branches of the arctan function 
in the inverse CAT ^ should be used (and a different to 7^ for the CAT). For each integer 
k, let us take T k = ((k — \)^, (k + \)^) and t k = fej. The solutions verifying conditions ^ 
at t k are u9°'(t) = (—\) k Ui{t) = U{(t — t k ) , i = 1,2. Define a pair of CAT and inverse CAT 

from K x T k into R 2 of the form: A {k) (x,t) = (jif^, JW^) = ( K > r ) and A lk)( K > T ) = (M)> 
where the fc-th branch of the arctan function has been used in ^4^- An unfolded version of 
the CAT, A, can be built by joining all the patches A^, defining an application that maps a 
complete trajectory x(t) of the harmonic oscillator into a trajectory k(t) of the free particle. A 
is periodic on t with period ^, although discontinuous. 

Other simple examples where this construction can be done are the damped particle and 
the damped harmonic oscillator, see [16] for details. 

3 The Quantum Arnold Transformation 

An arbitrary LSODE system ^ can be derived from the Lagrangian (we shall take A = for 
simplicity but the whole formalism can be developed with A ^ 0, see [3]): 

L = ^me f (x 2 - oj 2 x 2 ) , (9) 

and from this the Hamiltonian 

H = —e- f + 1 muj 2 x 2 e f (10) 

2m 2 v ' 
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is derived, which is known as the Generalized Caldirola-Kanai (GCK) Hamiltonian for a damped 
oscillator (see [3] and references therein). The case in which / = 7 and u are constants corre- 
sponds to the original Caldirola-Kanai Hamiltonian for a damped harmonic oscillator \17\ I18j. 
and whose corresponding Lagrangian was given for the first time by Bateman |19| , after a 
constraint to the Bateman dual system is imposed |20| . Canonical quantization of the GCK 
Hamiltonian leads to the time-dependent Schrodinger equation: 

The CAT A is a local (in time) diffeomorphism between the space of solutions of the LSODE 
system ^ and the space of solutions of the free particle. We would like to extend it to a unitary 
transformation A, the Quantum Arnold Transformation (QAT), between the Hilbert space of 



solutions (j)(x, t) of the time-dependent Schrodinger equation for the GCK oscillator ( 11 ) at time 
i, Tit, into the Hilbert space of solutions ip(K,r) of the time-dependent Schrodinger equation 
for the Galilean free particle 

ih^ = -—^ (12) 

dr 2m Ok 2 ' 
at time r, T~L^. The desired extension is given by: 

A: H t — > %r 

(f)(x, t) .— > <p(K, t)=A (<P{x, t)) (13) 



. i m 1 "2(*) ,.2 



A* y/u 2 {t)e 2 n ww U2 w x ^( X; f 



Here A* is the pullback of the CAT A, acting on functions (i.e. A*(f(x,t)) = r))). 
The QAT can be diagrammatically represented as: 



U G (r) 



u(t) (14) 



i 

where = T~Lq = Ti is the common Hilbert space of solutions of the Schrodinger equation for 
both systems at t = r = (we shall take, for simplicity, to = 0, as before), U(t) is the unitary 
time-evolution operator for the GCK oscillator and Ug( t ) is the corresponding one for the 
Galilean free particle. The map at the bottom of the diagram is the identity due to conditions 
([5]); otherwise a non-trivial unitary transformation appears (see [3]). 

From the commutative diagram, it is clear that A is unitary given the unitarity of the 
evolution operators. However, it can also be checked explicitly that the scalar product of two 
states in at a given time r is the same than that of the transformed states by A in T-L t at 
the corresponding time t: 
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/+oo 
dK <px(k,t)*<P2(k,t) 
-oo 



(i m 1 **2(*) 2 \ 
X|) e nff( t )^^ 1 ( Il ^) , (15) 
/ 

(/ ? m 1 "2( t ) T 2 N 

y / u 2 (t)e 2 n w ^ u 2 (t) x (j) 2 (x,t) 



X 
+oo 



dx (f>l(x,t)*(f)2(x,t) 



n,<P2)H t 



— oo 



where r, k, d« and the integration limits have been transformed according to the CAT. 

The QAT inherits from the CAT the local character in time, in the sense that it is valid only 
for t E T and tgT, although it can be defined for an arbitrary initial time to- To extend the 
QAT beyond T, we can proceed as in the classical case for the harmonic oscillator, considering 

the different branches of the inverse function of r(t), defining an unfolded QAT, A. 

It should be stressed that if in the different branches of the unfolded CAT proper solutions 
verifying §5§ are not used, changes in signs can appear, which result in changes in phases in the 
different branches of the unfolded QAT. This phenomenon is related to the Maslov correction 
(see for instance [21]). In fact, it can be checked that, for the case of the harmonic oscillator 
previously considered, A^(<f>(x,t)) = e lk 2 A((f>((— l) k x, t)). 

3.1 Symmetries of the quantum GCK oscillator 



By the results of Sec. 2.1 any classical LSODE system has the maximal number of Lie symme- 
tries, the ones of the free particle (SL(3, M)). The vector fields generating these symmetries can 
be computed by means of the CAT from the symmetries of the free particle. Even though the 
CAT is a local diffeomorphism, these vector fields are defined for all times (22] (the reason for 
this is that the only function appearing in the denominator of the vector fields is the Wronskian, 
which is never vanishing). 



The quantum symmetries of the free Schrodinger equation (12) are smaller than SL(3,M), 
only the Schrodinger group preserve^] ( fl2| [25] . The Schrodinger group is a semidirect product 
of the (centrally extended) Galilei group by the group SL(2, R) . It is generated by the basic 
operators: 

ifi d A d . 

K = K-\ T— , 7T = -Ifi— , (16) 

m on on 
in addition to their quadratic combinations (representing scale and non-relativistic conformal 
transformations) . 

The basic operators, together with their quadratic combinations, are constant of motion 
operatorg^J they are infinitesimal generators of symmetries of the free Schrodinger equation 
(12), and preserve the Hilbert space of solutions fll2h. 



We can import the symmetries acting on the free particle Hilbert space Ht (the Schrodin ger 



group) into the GCK Hilbert space Jit by means of the commutative diagram (14). The in- 
finitesimal generators of these symmetries are: 

P = -ihui- mx— , A =-xH ui— , (17) 

Ox W W m ox 

1 Symmetries of the Schrodinger equation are the quantum analogue of Noether symmetries, i.e. those pre- 
serving the Lagrangian. 

2 They verify that ^ + ^{H , ■] = 0, where H G = ~j^-§^2 is the Hamiltonian operator for the free Galilean 
particle. 
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and their quadratic combinations (see [3] for the complete list of operators). All these operators 



expand a representation of the Schrodinger group, are constant of motion operators for (11) 
and preserve the Hilbert space Tit- 

Despite of the local character of the QAT, and similarly to the classical case, these operators 
are well defined for all t € R, since the Wronskian W never vanishes. 

The Hamiltonian for the GKC harmonic oscillator can be expressed in terms of the basic, 
conserved operators: 

H = a(t)^ 2 + \mP{t)X 2 + 5{t)\{XP + PX) , (18) 
2m 2 2 

where the time-dependent functions a(t),/3(t) and 5(t) are given by: 

= «0-*^, m = "^">K (19) 

Again, since the Wronskian never vanishes and the basic operators X and P are also well- 
defined for all times, this construction remains valid for all times. 

Note that, since the basic operators X and P are conserved, H is not a conserved operator 
(unless uj = ujq is a constant and 7 = 0), and therefore an eigenvalue equation for it makes no 
sense (i.e., there is no time-independent Schrodinger equation). 

Instead, any linear combination of P 2 , X 2 and XP can be used to define an eigenvalue 
problem, and the corresponding eigenvectors would constitute a basis for the solution of the 
Generalized Caldirola-Kanai Schrodinger equation (see [3] and references therein). 

3.2 Applications of the QAT 



From the commutative diagram (|14|) and from ( 15 ) it is clear that the QAT is a unitary operator, 



and this has interesting and far-reaching consequences. Among them we can mention the 
possibility of importing operators (symmetries) from one system to the other, importing wave 
functions, scalar product, computing the time evolution operator, etc. We shall not discuss 
them here, referring the reader to O [T5l IM] for details. 

The QAT have many more and very interesting applications and generalizations. One of 
them is the application of the QAT to density matrices, transforming the Quantum Liouville 
equation of an arbitrary time-dependent quadratic Hamiltonian into the free one |15j . with the 
possibility of extending it and to Wigner functions and other quasi-probability distributions, or 
even the extension of the QAT to Lindblad type equations, in order to study dissipation and 
decoherence under the QAT point of view. It can also be applied to Bose-Einstein Condensates, 
to transform the time-dependent potential (oscillator traps with time-dependent frequencies) 
into a time-independent potential. 

The QAT can also be relevant in quantum inflationary Cosmological models, where quantum 
fluctuations of the scalar field are governed, under some assumptions, by a LSODE with a 
dissipative term. 

Among the possible generalizations of the QAT, we mention the relativistic case, as the 
quantum version of a generalized CAT for geodesic equations in a fixed background and with 
external forces, or as the quantum version of geodesic mappings, that transform geodesies of 
a metric into geodesies of a different metric (Beltrami Theorem). Other possible generaliza- 
tions refer to non- linear potentials (the Quantum Lie Transformation), to second order Riccati 
equations (see for instance [25]), or even to non-local potentials like in the Gross-Pitaevskii 
equation. 



7 



4 Extensions of the QAT 



In this section we shall introduce the main results of the paper: two immediate, although 
non-trivial, extensions of the QAT. The first one is to allow for more general Schrodinger 
equations than the GKC one. The second one is the transformation relating two LSODE- 
systems, denoted as the Arnold-Ermakov-Pinney transformation, obtained as the composition 
of the CAT corresponding to the first LSODE-system with the inverse of the CAT corresponding 
to the second LSODE-system, as well as the corresponding quantum version. 



4.1 Gauged Quantum Arnold Transformation 

The QAT, as it has been designed, transforms the solutions of the Schrodinger equation of the 
GKC oscillator into that of the free particle. It cannot transform, as it stands, the solutions of 
the Schrodinger equation associated with an arbitrary time-dependent quadratic Hamiltonian. 

The reason is the choice of solutions u\(t),U2{t) satisfying the canonicity condition ^ and 
the fact that the classical counterpart, the CAT, transforms the LSODE ^ into that of the 
free particle. 

Suppose that the quantum Hamiltonian is an arbitrary, time-dependent self-adjoint quadratic 
Hamiltonian: 

6 = " (4 + I) + k™®* ' (20) 

where T(t) and v{t) are real and > (to have a positive-definite kinetic energy). 

Then the transformation that maps solutions of the Schrodinger equation for this Hamilto- 
nian into that of the free particle is given by: 

(p(K, t) = G {(f>(x, t)) = A* (ju 2 (t) e~ 5 f ww ( 2$J ~ r (*)) * <f,( Xi t) \ , (21) 

where u\ (t) and ui (t) are canonical solutions (in the sense of j}) verifying the LSODE Q with 
/ = — - and uj 2 = fiv + T(^ — T) — T. The Wronskian of the two solutions is W(t) = nit). 

I I i m r(t) x 2 

Note that the transformation (21 ) is an ordinary QAT times an extra gauge term e 2 h ff W : 

i m r(«) 2 - 

G = A*(e* hW ® )A. (22) 

For this reason we shall denote this transformation Gauged QAT (GQAT). In the case where 
r(t) = 0, the GQAT turns into an ordinary QAT. 

Note that M(t) = -jjrk can be interpreted as a time-dependent mass, and therefore the 
Hamiltonian can be written as 

h 2 d 2 ( d IV 1 ,2 2 



H = -WW) 5? - mt)h W2) + 2 M («'> V • < 23 » 

where Q(t) 2 = fj,(t)v(t), and the classical damping coefficient and frequency are written now as 
/ = § and co 2 = Q 2 - T{§ + T) - f , respectively. 



4.2 The Arnold-Ermakov-Pinney transformation 

Although the CAT relates any LSODE system with the free particle, the dynamics of both 
systems can be very different]^] Thus, it could be interesting to relate directly two arbitrary 

3 Consider, for example, the harmonic oscillator with a bounded, periodic motion as compared with the free 
particle, with an unbounded non-periodic motion. 
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LSODE systems with similar behaviour, and this can be achieved by composing a CAT and an 
inverse CAT. 

More precisely, let A\ and A2 denote the CATs relating the LSODE-system 1 and LSODE- 
system 2 to the free particle, respectively, then E = AJ A2 relates LSODE-system 2 to LSODE- 
system 1. E can be written as: 

i?:lxT 2 ->8xTi 

(x 2 ,t 2 ) t-> (xi,h) =E(x 2 ,t 2 ). (24) 

The explicit form of the transformation can be easily computed by composing the two CATs, 
resulting in: 

*~m w ^-mh>- < 25 > 

where b(t 2 ) = L 2 satisfies the non-linear SODE: 

«2 (*i) 



Wo 2 1 

6 + ^ + ^ 6 =T#63 



2x h 2 W U 



(26) 



and where u\ refers to the i-th particular solution for system j; Wj, fj and ujj stand for the 
Wronskian and the LSODE coefficients for system j; and the dot means derivation with respect 
to the corresponding time function. 

For the particular case where LSODE-system 1 is a harmonic oscillator (ui(ti) = uq and 
fi = 0), this expression simplifies to: 

Wo 2 o 

b + f 2 b + U2b= -^-u 2 , (27) 

resulting in a generalization of the Ermakov-Pinney equation. For /2 = the Ermakov-Pinney 
equation (also known as Milne-Pinney) is recovered |26 [ 127 } [28]: 

b + oj 2 b=^, (28) 



representing a harmonic oscillator of frequency 102 with an extra inverse squared potential -§-. 
For ujq = 0, the Arnold-Ermakov-Pinney transformation reduces to the ordinary CAT, i.e. 
E = A. 

The Ermakov-Pinney equation is related to the Ermakov invariant |26t H] and appears in 
many branches of physics, such as Cosmology |29| . BEC [30J, etc. Its generalization to higher 
dimensions, know as Ermakov pairs (or system), appears in BEC \10\ 111) and what is known as 
Kepler-Ermakov systems [31j. 

The Ermakov-Pinney equation entails a kind of nonlinear superposition principle, in the 
sense that its solutions can be written in terms of the solutions yi(t),y2(t) of the corresponding 
linear equation (with loq = 0): 

b{tf = cm{tf + c 2 y 2 {t) 2 + 2c 3yi (t)y 2 (t), c iC 2-4=lo 2 . (29) 

The other way round, the general solution y{t) of the linear equation can be writen in terms 
of a paricular solution p(t) of the Ermakov-Pinney equation (|28[) as: 



y(t) = c lP (t)cos(co e(t) + C2), (30) 
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where c\,c% are arbitrary constantss and 8(t) = J p 2 dt'. Note that this equation is just (25) 



for W\ = W2 = 1, t\ = 8 fa), p = b, X2 = yfa) and x\ = yfa)/bfa) = c\ cos(u;oii + C2). As a 



result, the general solution of (28) can be determined from a particular solution p(t) using (30) 
and g9]|. 

The quantum version of the Arnold-Ermakov-Pinney transformation, E, can be obtained 
computing the composition of a QAT and an inverse QAT, to give: 

E: < } — < 

4>{x2,t 2 ) 1 — > (p(xx,ti) = E ((j)(x2, t 2 )) (31) 



(i m 1 ^{^2) 2 



The Quantum Arnold-Ermakov-Pinney transformation is a unitary map importing solutions 
of a GCK Schrodinger equation from solutions of a different, auxiliary GCK Schrodinger equa- 
tion which, in particular, can be the one corresponding to a harmonic oscillator. In that case 
the transformation is very similar to the one used in BEC, known as scaling transformation 
to transform the time-dependent potential (oscillator traps with time-dependent frequencies) 
into a time-independent harmonic oscillator potential [10^ lllj . Also, in that case (i.e. for 



/2 = 0, W2 = 1, /1 = 0, W\ = 1) equation (31) reduces to the transformation given by Hartley 
and Ray [32] (this was already given by Lewis and Riesenfeld in |4j). However, the Quantum 
Arnold-Ermakov-Pinney transformation allows to choose in a suitable way the auxiliary system 
from which the solutions may be imported. 
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